Noname manuscript No. 

(will be inserted by the editor) 



Microflares and the Statistics of X-ray Flares 

L G. Hannah 1 , H. S. Hudson 2 , M. Battaglia 1 , 
S. Christe 3 , J. Kasparova 4 , S. Krucker 2 , 
M. R. Kundu 5 , and A. Veronig 6 



< 



o 

(N 

£Q Abstract This review surveys the statistics of solar X-ray flares, emphasising the new views 

\ that RHESSI has given us of the weaker events (the microflares). The new data reveal that 

these microflares strongly resemble more energetic events in most respects; they occur solely 
within active regions and exhibit high-temperature/nonthermal emissions in approximately 
the same proportion as major events. We discuss the distributions of flare parameters (e.g., 
peak flux) and how these parameters correlate, for instance via the Neupert effect. We also 
highlight the systematic biases involved in intercomparing data representing many decades 
of event magnitude. The intermittency of the flare/microflare occurrence, both in space and 
\ in time, argues that these discrete events do not explain general coronal heating, either in 

^ ■ active regions or in the quiet Sun. 
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1 Introduction 

A solar flare is a rapid and transient release of energy in the solar corona, associated with 
electromagnetic radiation from radio waves to y-rays, local plasma heating to tens of MK, 
particles accelerated to GeV, violent mass plasma motions, and shock waves. In the largest 
events, a spectacular variety of phenomena can be studied in detail acr oss many wavelengt h 
ranges to try to understand the processes involved (see examples in iFletcher et al.ll201lh . 
The X-ray emission observed is of particular interest as it shows the accelerated particles 
and intense heating, and therefore provides the most direct insights into the physics of the 
basic energy release. Hard X-rays (HXRs), from about 10 keV to hundreds of keV , are pri- 
marily produced via thick-target bremsstrahlung (lBrownlll97ll : lKontar et al.ll2010h in which 
the coronal accelerated electrons are stopped instantaneously through Coulomb collisions 
with denser material in the lower solar atmosphere. Soft X-rays (SXRs), typically the com- 
ponent below tens of keV, are thermal emission (lines and continua) by plasma of a few to 
tens of MK. This emission results to some extent from heating at the site of energy release, 
but mainly fro m new coronal materi al evaporated from the chromosphere during the impul- 
sive phase (see IFletcher et al]|201 ll) . Given the wealth of X-ray flare data, it is possible to 
study large samples of these events, allowing the statistics of the events to provide clues to 
the underlying processes behind the emission. This article reviews such flare studies, the re- 
sults, interpretation and limitations, from a predominantly X-ray viewpoint. In particular, we 
discuss the advances that have been made in extending this analysis to weaker HXR events 
(microflares) us ing the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) 
I Lin et al.ll2002h . 

One of the most remarkable features of flare surveys is the range of magnitudes these 
transient phenomena span. The observed SXR flux extends over five orders of magnitude 
(see Figure II. Il l and over nine orders of magnitude in energy (see Figure 12.11 we return 
to discuss these Figures in detail in Section [3). The observations come from many space- 
craft: e.g., Orbiting Solar Observatory (OSO-3) and (OSO-7), the Solar Maximum Mission 
(SMM), the International Cometary Explorer (ICE, a.k.a. ISEE-3), the Compton Gamma- 
Ray Observatory (CGRO), Yohkoh, the Solar Heliophysical Observatory (SOHO), GRANATI 
WATCH, and the Transition Region and Coronal Dynamics Explorer (TRACE). The dis- 
tributions of these and other flare parameters, such as emission at other wavelengths and for 
other durations, have the further striking property that they all can be well represented by a 
power law of the form 

f(x;a) = Cx~ a , (1.1) 

where / is the probability density function (PDF; this is often called the flare frequency dis- 
tribution) of the flare parameter x, a > is the power-law index and C is a scaling constant. 
In the case of the flare's energy U, the quantity f(U)dU is the fraction of events per unit 
time releasing energy between U and U +dU. This power-law natur e of the frequenc y dis- 
tribution has long been observed, first noticed in solar radio bursts bv lAkabanel d 19561) . The 
fact that so many of the flare characteristics have s uch a distribution is th ought to arise from 
the corona being in a self-organized critical state dLu&Hamiltor]|l99ll) . discussed further 
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Fig. 1.1 Frequency distribu- 
tion of the peak 1-8 A SXR 
GOES flux of 49,409 flares 
between 1976 and 2000 
iVeronig et al]|2002al) . Note 
that this study did not in- 
clude the A-class events 
(<1(T 7 W m" 2 ). In addi- 
tion, there was no attempt 
at background subtraction, 
so the resulting value of a 
should be regarded as an up- 
per limit (background pre- 
dominantly affects smaller 
flares, artificially steepening 
the distribution). 



in Section 13.31 Here, the underlying energy release from the largest to the smallest flares 
occurs as an "avalanc he," similar to the theory for earthquake occurrence or stress relief 
in general dKatzll 19861) . These properties have brought forth interesting interpretations. The 
power-law distribution means that small flares are more numerous. 

Could these weaker events at the currently unobservable "nanoflare" magnitude level 
(about 10~ 9 times the energy of a large flares) directly constitute the h ot corona, thu s ex- 
plaining its heating? This is often called Parker's nanoflare mechanism JParkeil 1 1 9 8 8h . The 
answer to this question requires knowledge of the total power that is contained in the flare 
distribution, i.e. 

P = / f{U; a)UdU = - [U m J- a - t/ min 2 - a . (1.2) 

Since U m i n <S (/ max by definition, we can see that in the case a > 2 in Equation dTTSb the 
low-energy half of the distribution would contain the most energy dHudsonlll991bl) . In this 
situation the smallest flares fit the requirement for heating the corona as they have a high 
occurrence rate and release more net energy than the large flares. Therefore, to understand 
the role of flares in coronal heating, the power-law index a has to be accurately estimated; 
this is problematic as the energy is inferred from observations and thus subject to potentially 
large errors and biases. In addition, Equation i ll. 2b assumes that the distributions continue 
into the unobservable low-energy range. Even if the deduced a were large enough (i.e., 
larger than two), that would not automatically mean that small flares heat the corona. For 
further proof we would require observations showing the power-law index to maintain its 
value down to the distribution's physical limit, rather than its observational limit. 

Between the largest flares and these theoretical nanoflares lie the microflares (nomi- 
nally with energies of order 10~ 6 times those of large flares), and it is with these events 
that RHESSPs HXR observations have made major advances. RHESSPs continuous spec- 
tral coverage, down to 3 keV, uniquely allows the study of the transition between the thermal 
and nonthermal emissions with the same instrument. In particular, RHESSPs view of their 
nonthermal characteristics has allowed us to study smaller-scale active region events with 
simpler structures than the major events that often attract the most attention. Perhaps these 
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can help us to isolate the essential physics, which could then be applied in more complex 
situations. 

We begin this article with a general overview of flares and flare-like brightenings across 
the range of magnitudes (Section 0. This discloses similar properties between major and 
minor events, but systematic differences do occur. In Section[3] we learn what we can from 
these differences about the physics, taking advantage of the very large numbers of events 
in major surveys, including RHESSI's new HXR views. We also discuss in this section the 
biases that arise in the statistical surveys and possible methods for obtaining the unbiased 
intrinsic distributions. In Section [3~3l we briefly discuss how the power-law nature of the 
flare parameters arises. Conclusions and discussions are given in Section[4] 



2 From major to minor flares 

2. 1 Flare classification & general properties 

The most powerful ordinary flares have energies estimated at above 10 33 ergs and present 
a spectacular range of phenomena, easily observed across the wavelengths. The first flare 
observed was a powerful event i n 1859, detectable through its small, i ntense white-ligh t 
emission patches as described by ICarringtonl ll 1859h and corroborated bv lHodasonl d 18591) . 
Remarkable terrestrial effects accompanied this flare and also followed it after an interval of 
half a day. This event anticipated much of the complexity of flares as we know them today, 
but it was no t until the 1940s that "flare" was a ccepted as the term to describe these transient 
phenomena dNewtonlll943l : lRichardson|ll944h . Events with total energy ab out a millionth 
smaller than large fla res (about 10 27 erg), became known as "microflares" dSchadee et al.l 
1 19831 ; I Lin et al Parker hypothesized that even smaller flares, "nanoflares," with en- 

ergies of order one billionth of lar ge flares or a bout 10 24 erg, could be the basic unit of a 
localized impulsive energy release jParkeilfl988h . 

Quantitative flare classification is based on the 1-8 A SXR flux observed by GOES. 
Large flares have Xn-class, indicating a peak flux of n x 10~ 4 Wm~ 2 , the largest events 
being above X10. This classification decreases through the decades of M, C and B-class 
flares down to the smallest An-class events with n x 10~ 8 W m~ 2 and the sensitivity limit of 
the detector. The classification of flares and the associated range of GOES fluxes is shown 
in Table ITTI The largest GOES flare was SOL2003-11-04T19:53 (X17.4), which saturated 
the detectors at 18 x 10 -4 Wm~ 2 (i.e., class X18). It is es timated that this flare was an 
X28 event with range X25 to X31 ( iKiplinger & Garciall2004l) . In terms of GOES classifica- 
tion, the Carrington fla re appears to have been a large soft X-ray flare of magnitudfe >X10 



(Tsurutani et al. 2003) and like many large flares resulted in a major geomagnetic storm 
Cliver & Svalgaardll2004l) . This region also produced one of the larg est solar energetic pro- 



(Liiver & svalgaard z(JU4 ). inis region also produced one ot me largest solar energetic pro- 
ton fluences at the Earth in the last 500 years dMcCracken et al.l200ll) . These most energetic 
flares occur in active regions, often when new flux emerges into an already-complex mag- 
netic structure. From the start of the GOES observations in 1975, to the start of cycle 24 in 
2009, 22 flares >X10 had been observed. In total, 359 GOES X-class flares, 4708 M-class 
32784 C-class and 11558 B-class flares occurred between 1976 and 2000 jVeronig et al.l 
l2002ah . Although more small flares are expected due to the inherent power-law nature of the 
flare distribution, they are hidden by the higher background during active times. The major- 
ity of large flares seem to occur during the peak of the solar cycle headin g into the deca y 
phase, though this behavior can vary dramatically from cycle to cycle (e.g., lHudsonll2007l) . 
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Fig. 2.1 The energy dis- 
tributions for solar flares. 
The nonthermal energy 
distribution is shown for 
large flares >25 keV ob- 
served with S MM/H XRBS 
ICrosbv et all Il993h . mi- 
croflares >8 k eV from 
CGRO/BATSE iLin et al l 
I2001ah and microflares 
> Eq (above the low en- 
ergy cutoff) from R HESSI 
iHannah et al] I2008al) . The 
thermal energy distribution 
is shown fo r microflares 
with R HESSI iHannah et all 
l2008al) and Yoh koh/SXT 
iShimizj 119951) and 
EUV nanofiares with 



TRAC E iParnell & Jupt 



200C ; lAschwanden et al 



20001) and S O HOIEY I 
)Benz & Kruckeri 120021) . 
This figure is deceptive 
as it is comparing energy 
distributions of different 
flare energy components, 
each involving different 
instrument and selection 
effects, and were obtained 
over different periods of 
different solar cycles. 



Table 2.1 Different flare classifications and the associated ranges of GOES flux (SXRs) and HXRs. 
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"Active Region Transient Brightening" observed in SXRs jShimizull 19951) 



These two features can be seen in the GOES flaring rate shown in Figure [2~4l discussed in 
Section lZZTI 

Flare phenomena encompass all of the accessible wavelength ranges of electromagnetic 
radiation, a s well as the emissio n of neutral particles, gaseous ejecta, and large-scale shock 
waves (e.g. lFletcher et al.l201 if) . In the most energetic flares, all of these phenomena appear, 
considerably more intensely than in smaller flares and re gardless of the detailed physics; 
this is the so-called big-flare syndrome (e.g.. iKahlerll 19821) . This makes these events help- 
ful diagnostically. This property suggests that if a particular phenomenon (e.g., white-light 
continuum) is not observed in a weaker flare (as not distinguishable from the background), 
the process may still function in the same way. Non-scalable properties might exist; one 
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Fig. 2.2 The association rate of CMEs with 
flares as a function of the flare peak GOES flux 
lYashiro et aT] I20O6I) . Only flares from C-class 
(KT 6 W m~ 2 ) to X-class (>1(T 4 W nT 2 ) are in- 
cluded in this survey. Reproduced by permission 
of the AAS. 



that is often cited is the maximum energy attained by accelerated particles. An emission 
such as ^°-decay y -ray emission might thus have a non-linear threshold dependence; see 
IVilmer et al] J201of) . Of course, given that there weak flares are considerably more numer- 
ous, exceptional cases can ofte n be found, such as m icroflares with nonthermal emission to 
remarkably high energies (e.g. jHannah et al]|2008bl) . 

2.1.1 Flare relationship to CMEs and SEPs 

There has been a long debate whether ejection of material from the corona, namely a Coro- 
nal Mass Ejection (CME), triggers a flare, or vice-versa, but it is now generally accepted 
that they are both consequences of coronal energy release through the reconfiguration of 
the magnetic field. The issu e of whether a flare o r CME is produced will depend on the 
magnet ic field configuratio n dWang&Zhandl2007l) . Using a large survey of nearly 7000 
CMEs, lYashiro et alTj2004l) found that all the most energetic (>X2) flares have CMEs; see 



Figure [2T2l ( Yashiro et al The largest and fastest CMEs are generally associated with 



such large flares. The rate for C-class flares is about 20% and for M-class about 40%. Of 
course, there may be a sensitivity issue which limits the ability to observe the small CMEs 
associated with microflares. Studies of flares that do not have associated CMEs have foun d 
these events to have a correlation between peak intensity and flare duration jKav et al.l2003h . 
No such correlation was found for CME flare events, suggesting a physical difference in the 
time development of a flare when a CME is involved. CMEs appear to have broadly the same 
range of characteristics whether or not they are associated with a fl are, but the fastest and 
broadest CMEs are always associated with the most energetic flares jVrsnak et alj|2005h . In 
terms of the energy partition between fl ares and CMEs, the CME kinetic energy ranges may 
be comparable to the total flare energy femslie et al.ll2005h . 

The ejection of highly accelerated (up to several MeV) electrons, protons and heavy ions 
from the Sun into interplanetary space, called Solar Energetic Particles (SEPs), are some- 
times associated with flares (e.g. lReamesll 19991) . Impulsive SEP events demonstrate a rapid 
burst of electron enhancement often related to the impulsive HXR phase of flares, suggesting 
that the same coronal energy release is responsible for the flare and SEPs. These energetic 
electrons therefore provide an important diagnostic of flare particle acceleration, especially 
when combined with HXR flare observations. Type III radio bursts are also produced by en- 
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Fig. 2 .3 (Left) The locatio ns of 24,097 microflares observed with RHESSI between March 2002 and March 
2007 iChriste et all2008dl All of the confirmed events could be associated with active regions. Reproduced 
by permission of th e AAS. (Right) The locations of one day of RHESSI microflares from active region 10456 
I Stoiser et al. 2007). Of 53 events found, 24 could be imaged. The centroid positions (3-6 keV) are indicated 
by crosses and overplotted on the MDI (Michelson Doppler Imager of SOHO) magnetogram. 

ergetic electrons escaping from the Sun into interplanetary space, see Section l2.2.3l Gradual 
SEP events show a slower enhancement of energetic protons, thought to be accelerated by 
a CME-associated shock in interplanetary space rather than in the coronal en ergy release. 
Although these proton events are not directly related to flares, iKiplingeri Jl9951) showed that 
flares that were associated with 10 MeV proton events predominantly demonstrated pro- 
gressive hardening of the HXR photon spectra into the decay phase of the HXR emission. 
In contrast, the other population of HXR events not associated with SEPs exhibit softening 
spectra in the decay ph ase. A similar result w as found for X-class flares during the January 
2005 solar storm event dSaldanha et"a7]|2008l) . In these flares, four out five showed spectral 
hardening with RHESSI and were associated with interplanetary proton events. A study of 
37 events with RHESSI found that the majority of events that demonstrated HXR spectral 
flattening (12 out of 18) produced SEPs (detected with GOES >10 MeV proton data and 
WIND/3DP 0.1-1 MeV proton and 30 - 500 keV el ectron data) and all without flattening 
(19 flares) did not produce SEPs dGravson et aDl2009h . 



2.2 Microflares 

2.2.7 Association with active regions 

"Subflares" had always been known to the Ha observers (e.g., ISmith & Smith! 1 19631) as 
events that were small in area but not nece ssarily faint. The t erm "microflare" was intro- 
duced in the 1980s bv lSchadee eTail dl983h for SXRs, and bv lLin et alJ Jl984h for HXRs. 
The presence of X-rays, the occurrence of subflares wi th "brilliant" Ha classification (high 
inten sity), and the discovery of faint microwave bursts dGopalswamv et ail 1 9941 ; rGarv et alJ 
1 19971) in association with SXR microflares made it clear that the basic flare physics extended 
over a wide magnitude range. Flares generally (including subflares and microflares) show 
both thermal and nonthermal emission, indicating the presence of particle acceleration and 
plasma heating. In terms of GOES classification these events are typically A- and B-class 
events, down to and beyond the sensitivity limit of the GOES detectors. 
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Fig. 2.4 (Left) Occurrence rate of RHESSI microflares (<C1 events) {Christe et al. 2008a), showing a clear 
decline over the waning years of Cycle 23. Reproduced by permission of the AAS. (Right) The number of 
GOES flares per year per flare class for X-class (top) through M, C to B-class (bottom). The microflares seem 
to be anti-correlated with the larger flares, which would be expected from the increased back ground duri ng 
times of high activity. Note that B-class flare data was only available from 1984 (data from Milligan 2009). 



The HXR observations prior to RHESSI normally suffered from poor sensitivity due to 
lack of access to lower hard X-ray energies (<15 keV), hence these observations missed 
most of the microflares. The problem was the use of thick fixed attenuators to reduce the 
excessive low energy counts in large flares. One exception to this was the HXISQ imager on 
the Solar Maximum Mission (SMM) dSchadee et alJll983h . which covered the 3 — 30 keV 
range with a restricted field of view and very small detectors, thus avoiding saturation. Be- 
cause of this, the main HXIS results were in a spectral region comparable with the GOES 
data. At lower energies (0.25 — 4 ke V), the Soft X-r ay Telescope (SXT) on Yohkoh found 
microflares to occur in active regions dShimizu|[T995l) . . 

The first observation of HXR microflares was with balloon-borne detectors that could 
observe only down to about 15 keV due to atmospheric absorpti on, and which had high 
sensitivity from the use of large-area, low-background detectors dlin et ai] ll984» . Subse- 
quently with Yohkoh/HXT, emi ssion at 14 — 23 keV was found to be associated with SXR 
brightenings in activ e regions i Nitt3 1997t) . Microflares were also detected at 8 — 13 keV 



with CGRO/BATSE dlin et al.ll2001bh . However, it was not until RHESSI that the detailed 
associations of these HXR microflares could be readily investigated. 

RHESSFs introduction of a movable attenuating shutter system dlin et alj2002h allowed 
the attenuators only to be deployed when the detectors were saturated. With no shutters, the 
detectors become saturated once solar emission (either background emission from active re- 
gions or flares) is approximately greater than GOES CI level. RHESSI microflares are there- 
fore flares that occur when the attenuating shutters are out; they are typically sub-C-class 
events, < 10~ 6 W m~ 2 in the 1 — 8 A GOES band. The full view and sensitivity of RHESSI is 
available to such events, providing the first comprehensive HXR view of microflares, with 
both imaging and spectroscopy. In addition to this, RHESSI introduced Ge detectors with 
resolution much superior to the scintillation counters previously used for HXR observations, 
allowing the properties of steep spectra to be accurately determined. 

RHESSI microflares appear to be exclusively localized to active regions. Figure [2~3l (left 
panel) shows this explicitly with the positions of 24,097 events, all of which are identified 



Hard X-ray Imaging Spectrometer. 
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as coming from active regions dChriste et al.ll2008ah . These events were identified from an- 
alyzing all count-rate spikes in the RHESSI 6—12 keV energy range between March 2002 
and March 2007, at times when the attenuating shutters were out. A total of 25,705 events 
were found, though only 24,097 events could be imaged to give positions on the solar disk. 
The events not counted include some of the faintest true microflares, but are predominantly 
misidentified particle events (which gave a non-solar disk position) or events close to the 
rotation axis of RHESSI (events close to where the spacecraft rotation axis is pointing will 
have no modulation in their detected signal and hence no spatial information can be recov- 
ered). There is repeated microflaring from active regions, thus tracing out the location of 
a region as it moves across the disk. The microflaring rate varies greatly between active 
regions. Detailed RHESSI imaging of all the microflares from one active region compared 
to MDI magnetograms show s that the microflar e emission occurred at multiple locations 
throughout the active region dStoiser et alj|2007l) . though many originated in one localized 
area of the region (Figure [231 right panel). Another sampl e of microflares showed that they 
mostly occur near magnetic neutral lines in active regions dliu et al]|2004l) . 

As expected, then, Figure 12.41 (left panel) shows that these events tend to follow the 
solar cycle in their occurrence frequency. We note that it is difficult to follow the microflare 
rate over a whole solar cycle because the solar X-ray background level also follows the solar 
cycle. Fewer microflares are observed during periods of high solar activity as they are hidden 
by the high background from active regions. This can be clearly seen in the number of GOES 
B-class flares (right panel of Figure l2.4b . Their rate appears unexpectedly to be almost anti- 
correlated with the solar cycle. This observational bias means that the microflaring rate is 
only accurately known during times of lower activity/solar background. 

As the HXR microflares observed by RHESSI are inherently active-region phenom ena, 
it seems entirely reasonable to associate them with microflares dSchadee et alJ 19831 ) and 
the active-region transient brig htenings (ARTBs) seen in soft X-rays lshimizulll995 ). This 
had been p reviously po stulated dTsuneta & Lemenll 19931) and early evidence suggested it to 
be correct JlSlittall 1 997h . but the confirmation was only possible with RHESSI's hard X-ray 
sensitivity. 



2.2.2 Physical properties 

The X-ray time profiles of typical microflares (examples given in Figure 12.5b demonstrate 
a similar structure to those in typical large flares: a short impulsive burst at higher energies 
(HXRs) followed by a more gradual thermal emission at lower energies (SXRs). Microflares 
also demonstrate other temporal phenomena seen in large flares like more gradual emission 
or pre-impulsive behavior though the timescales are generally shorter (seconds to minutes 
instead of several to tens of minutes). 

The first HXR spectra of microflares (top panel of Figure [2T6t showed only a power-law 
component, character istic of nontherm al emission, due to the detector being sensitive only 
above about 13 keV dlin et aljfl984h . RHESSI spectra (bottom row of Figure 12.6b can be 
obtained with moderate energy resolution down to a few keV. This allows the thermal com- 
ponent, and the transition from it to the nonthermal component, to be seen. These spectra 
(from the time of peak emission just after the impulsive phase) are similar to those seen in 
more energetic flares, except that (a) the thermal emission tends to be at lower temperatures 
and (b) the nonthermal component tends to have a steeper spectrum (fewer electrons accel- 
erated to high energies). In the bottom left panel of Figure [2761 we have a t ypical RHESSI 
microflare with a fitted temperature of about 13 MK. The Fe feature (e.g., |phillipsll200i) 
is clearly in evidence at about 6-7 keV and a nonthermal spectrum above 8 keV with an 
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Fig. 2.5 The GOES 1 — 8 A lightcurve (top panels) and RHESSI lightcurves for some example microflares. 
(Left:) A single B-class microflare (Bl with pre-flare background subtracted) which shows the classic flare 
time profile: impulsive at higher energies with a slower rise and more gradual fall at lower energies. This 
event also shows some pre-flare emission. (Right:) two A-class events (Al and A3, background-subtracted) 
both of which have no emission above 1 2 keV observed by RHESSI; their nonthermal emission hidden by 
the background. Note that the background level in these microflares is predominantly solar for GOES, but 
terrestrial and instrumental in origin for RHESSI. 



g-5.15 power law. The other RHESSI spectrum shown (bottom right panel of Figure [276b is 
of an anomalously hard microflare, which despite being a GOES A2 event (near the limit of 
the GOES sensitivity) demonstrates st rong nonthermal emi ssion to relatively high energies, 
just as in much more powerful events I Hannah et al.l2008bh . One should note, however, that 
the early microflare HXR spectra of iLin et alJ 1 1984 ) were more sensitive above about 30 
keV since their shielded detectors were able to reduce the background clearly evident in 
the RHESSI spectra. These Lin et al. (1984) microflares also demonstrate steep nonther- 
mal emission of /(e) °c with y = 4 to 6 that continue to energies normally obscured by 
background in RHESSI observations. 

RHESSI microflares are not spatially small and typically show an elongated loop-like 
struc ture, with higher-energy HXR emission from the footpoints at the ends of the loop 
(e.g jKrucker et alj|2002l ; iLiu et al.|[2004l ; IStoiser et alhocm lHannah et al.ll2008bh . Figure 
!2.7l shows an example of this overlaid on an EUV image from TRACE. We interpret the EUV 
loop as the cooling post-flare loop which was first seen at higher temperatures in the RHESSI 
image. This hot material could be driven (evaporated) from the chromosp here by the energy 
deposition of the nonthermal electrons penetrating to the loop footpoints dFisher et al.l 19851 : 
lAbbett & Hawlevll 19991 : [Brosius & Phillipsll2004ri and producing HXR sources there. Sev- 
eral RHESSI microflares ha ve also been observed with the HXR loop connecting two bright 
Ha kernels dLiu et al.l200i) , another possible signature of the energy deposition of nonther- 
mal electrons in the lower atmosphe re. Previous work ha d shown the association between 
microflare SXR and Ha emissions ( Shimizu et al. I l2002l) . consistent with the relationship 
known from ordinary flares ( iThomas & Teskdll97ll) . 

Fi gure [2~8l show s the fitted lengths and widths of the thermal emission (4-8 keV) for 16 
seconds about the tim e of peak emission in 6-12 keV for almost 19,000 RHESSI microflares 
dHannah et a71l2008d) . These parameters are from simple forward-fitted image models and 
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Fig. 2.6 (Top four panels:) the first HXR microfiare spectra from a balloon-borne instrument (L ineLalJ 
1984), and (bottom two panels: spectra from RHESSI showing a typical m icrofiare ^Hannah et afl l2008a) 
(left) and one with a particularly hard spectrum (right) I Hannah et al. 2008b). The red and blue lines are the 
thermal and nonthermal components of the model, respectively. In the RHESSI spectra, the upper limit of 
the nonthermal component is hidden by the instrumental background (the large errors bars in the bottom-left 
spectrum and the grey histogram in the bottom-right spectrum). Reproduced by permission of the AAS. 
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Fig. 2.7 Overlay of RHESSI ther- 
mal (red) and nonthermal (blue) 
contours on TRACE, showing the 
typical flare and microfiare sce- 
nario of (left) footpoints deposit- 
ing their energy to heat chromo- 
spheric material that expands up- 
wards, filling the loops connecting 
them. (Right:) his material then 
cools, emitting brightly in EUV. 
Updated vers i on of figure from 
lKruckeretai]<2002l) . 



show that the microflares have a clear tendency to be elongated, as we wo uld expect from 
their identification with the loops that characterize SXR microflares (e.g., IShimizu|[l993) . 
The distribution of the loop length has a peak well above RHESSPs angular resolution, so 
we conclude that these microflares may not be inherently small-scale. Many ordinary flares, 
though more energetic, have dimensions in this range. In addition, these l engths are not 
correla ted with the microfiare peak GOES flux or peak RHESSI emission dHannah et alj 
l2008al) . 
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Fig. 2.8 Loop arc length and width distribu- 
tions for 18,656 microflares imaged by RHESSI 
indicating that they are predom inantly elongated 
structures tHannah et aUI2008al) . Thus microflares 
are not necessarily spatially small. Reproduced by 
permission of the AAS. 



2.2.3 Associated ejecta and radio emission 

Radio signatures generally guide us to the corona, since the opacity (and the minimum 
height of fo rmation of the em ission) increases in height almost monotonically with wave- 
length (e.g. Jwhite et aLlfeoiOh . Millimeter-centimeter wavelengths (microwave GHz emis- 
sion) originate mainly in the lower atmosphere, meter-waves (100s MHz to MHz radio 
emission) in the lower corona, and the longest wavelengths right out into the solar wind. 
The longer-wavelength radio observations th us often show ejecta, especially via nonthermal 
plasma signatures (e.g., IPick & Vilme3l2008h . 

The microflare gyrosynchrotron emi ssion is closely correlated in time to the HXR emis- 
sion dOiu et al]|2004l ; lKundu et alj|200rj) , indicating that they are both signatures of the ac- 
celerated flare electrons. The previous observations of SXR microflares als o demonstrated a 
correlation with microwave emission dGarv et aHll997l : [Nindos et al.ll999h . The microwave 
spe ctral index for m icroflares is found to be flatter (harder) than the RHESSI HXR spec- 
tra dOiu et af]|2004) . This is thought to be due to the microwa ve emission originating from 
electrons with higher energies than those producing X-rays dNitta & Kosugil 19861) . This 



discrepancy in spectral indices is also observed in larger flares dSilva et alJl200Ch though 



it is possibly greater in microflares. This may be due to the HXR spectra being typically 
steeper in microflares than large flares. A B-class microflare was observed to have nonther- 
mal microwave em ission from MeV-energy electrons but showed no HXRs at lower energies 
dRaulin et al.lll999T) . This behavior could result from trapping of the lower-energy electrons 
during the peak of the microwave emission from the higher-energy electrons. 

Imaged microwave emission using the Nobeyama radioheliograph (NoRH) for RHESSI 
microflar es indicated similar ly sized and separated footpoints in both HXR and microwave 
emission dKundu et al]|2006h : see Figure [Z9I The higher-frequen cy microwaves come from 
the footpoints, and the lower-frequency from the connecting loop dKundu et al .120051 1 20061) . 

RHESSI microflares are often associated with Type III radio bursts, produced by elec- 
tron beams escaping from the corona and often subsequently detected as SEPs. These are 
the " fast drift" bursts identified with weakly relativistic electron beams (e.g., IWild et alJ 
Il963l) . The HXR emission in these cases still seems to be a signature of accelerated elec- 
trons reaching the lower atmosphere at loop footpoints, although the radio emission comes 
from accelerated electrons of similar energies moving outwards from a postulated coronal 
acceleration site. An example of six microflares identified with RHESSI and their associated 
MHz radio emission observed with WEVD/WAVES is shown in Figure IXTOl Here we can see 
the time correlation of the HXR emission and the high-frequency radio emission. The mi- 
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Fig. 2.9 Five microflares (one per column) imaged in different HXR energ y bands (increasing from top row 
to bottom) with RHESSI (contours), overplotted on NoRH 17 GHz images iKundu et al.ll2006h . The bottom 
microwave image is background-subtracted. Similar-sized footpoints and loop structures are seen in both 
HXRs and microwaves. 



croflare with the largest HXR flux produced the brightest radio burst. The association with 
Type III bursts establishes that electrons accelerated in the microflares have access to open 
field lines (or field lines that extend into the upper corona), not requiring an eruption prior to 
the event . An associat i on be tween Type III bursts and small SXR flares had previously been 
found bv lFuerst et all Jl982h . 

The presence of open field lines in microflare s has also been seen with the association 
of jets of material flowing out of the flare region dShimoio et al. Il 19961). For RHESS I HXR 
microflares, such jets have bee n observed in EUV with TRACE dChriste et al ]|2008bh and in 
EUV and SXRs with Hinode dChifor et al.ll2008h . In the former case, multiple microflares 
and type III burts were obse rved, occurring ever y few minutes, with EUV jets associated 
with the largest microflares dChriste et alJl2008bh . In the latter case, the recurring jets oc- 
curred on a timescale of hours and were attr ibuted to chromosp heric evaporation flows due 
to recurring coronal magnetic reconnection dChifor et alj2008l) . 

This association with jets is also well observed in soft X-rays dStrong et al.ll992l : IShibata et al.l 

Il992h . Indeed all such jets, including those in the quiet Sun, appear to have loop brightenings 
at their base. A jet appears as a collimated flow, implying the prior existence of large-scale or 
even open magnetic fields. Type III radio bursts, known to be produced by electron beams 
escaping from t he corona into the solar wind, have a strong association with the soft X- 
ray jets as well dAurass et alJI 19941 : IKundu et al.ll 1995h . Thus the microflare/jet events also 
have thermal and nonthermal attributes just as major flares do, but the jet is relatively more 
prominent in the microflare domain. 



14 



Hannah et al, 



i ; - 








lit /^V\A/1 hhf~*\tAfrP\* / ■ ■ An 
Mi uV 


2 H 






- 


Si li 




'i^Jlii ■■ ' ■■-=■ " — ■ - -■ .-...»—. , 


4 1 




5 i 








6 1 




"I 




7 - 



Fig. 2.10 The lightcurves and radio spec- 
trograms of six R HESSI microfl a res from 
u _ 15 keU July 19, 2002 IChriste et alj l2008bl) . 

The top two panels show the GOES 
lightcurves in 1 — 8 A (with a B8 back- 
4 - 7 kev ground level) and 0.5 — 4 A. The next 
three panels show the RHESSI spec- 
trogram with night-time background re- 
trace(i95a) moved and lightcurves in 12 — 15 keV 
(nonthermal) and 4 — 7 keV (thermal). 
" ,, °"' 5 The EUV lightcurve from TRACE and 
a radio spectrogram from W7AD/WAVES 
are shown in the bottom two panels. The 
radio emission starts at around 14 MHz 
and drifts down to 2 MHz. It is well cor- 
nier 2 ' related with the RHESSI HXR lightcurve. 
Reproduced by permission of the AAS. 



2.3 Nanoflares and non-active-region phenomena 

We have seen above that the RHESSI microflares occur only in active region s and have a 
strong kinship with "o rdinary" flares. How do these r elate to X-ray bright po ints ( iGolub et alj 
119741) . network flares jKruckeret aljl997h . plumes dAhmad & Webblll978b . or any of many 



more types of weak flare-like brightenings seen both inside and outside the active regions? 
Note that EUV observations such as those of SOHO/EYT or TRACE may be much more 
sensitive than SXR (o r HXR) observations of transient features at ordinary coronal temper- 
atures of order 1 MK jPorter et alJl987ljKrucker & Benzlll99ot) . Many of these phenomena 
have been related to Parker's nanoflares and to the problem of coronal heating. The nanoflare 
description of these EUV events is justified, in the sense that an EUV brightening may have 
a much smaller event energy than a microflare observed in SXRs, HXRs, or microwaves. 
However it is often difficult to place these different sorts of observations on a uniform en- 
ergy scale. We deal with this at length in Section [3~T1 below. 

Outside of active regions, the most prevalent X-ray features of SXR images are the X- 
ray bright points (XBPs), originally discovered with rocket-bome X-ray image rs in the late 
1960 s and studied statistically with data from the X-ray telescope on Skylab dGolub et alj 

Eazl). 

Subsequent observations had sufficient resolution to de tect hot (about 2 MK) loop-like 
(abou t 10" in length) prop erties and footpoints in the XBPs dMoses et alJl994l ; IStrong et alj 
ll992l : lKotoku et alJl2007T) . pointing to a possible relationship with flar e physics. Typicall y 
hundreds of XBPs are visible, uniformly spread across the solar disk dGolub et alJll974f) . 
their number changi ng little once observational bias, d ue to the dominant active region 
emission, is removed l lHara & Nakakubo-Morimotoll2003l) . This spatially-uniform and time- 
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Fig. 2.11 Observations and upper limits of the quiet-Sun hard X-ray flux found with RHESSI. (Left:) the 
RHE SSI upper limi ts to the quiet Sun emission at 3 — 200 ke V in comparison to p revious observations 
(Hannah et al. 2010), including the first RHESSI quiet Sun limits (Hannah et al. 2007b). (Right.) comparison 
of the 3 — 6 keV flux observed with RHESSI to the 1 — 8 A GOES emission for the quiet Sun and offpointed 
micro flares. The Ymp and y ps indicate the index of power-law fits to the data, the dashed and dashed-dotted 
lines iHannah et al.l2 007b). Reproduced by permission of the AAS. 



invariant distribution across the solar surface is expected as the XBPs are associated with 
the magnetic fields related to surface convective flows (the granulation and supergranulation 
of the quiet photosphere) rather than with the active-region fields responsible for sunspot 
fields. X-ray bright points are associated with connectivity changes in the network magnetic 
field. The underlying fundamental physical processes could be the same as for active-region 
microflares. Often, XBPs may persist for hours to days, unlike the transient active region 
flares, and thus their beh avior may resemble that of quiescent active-region loops (e.g., 
lYoshida & Tsunetalll996h . In the absence of solar active regions, the XBPs and the corona 
of the quiet Sun provide a basal X- ray emission level. Jets are sometimes observed from 
XBPs , particularly in polar regions dShibata et alJI 19921 : IShimoio et aDl2007l : ICirtain et alj 
l2007r .l ISavcheva et alj |2007) where there is easy access to open field lines. 

There have been no HXR observations of XBPs. Such emission therefore must be below 
the sensitivity limits of current spacecraft, which have been optimized to observe brighter 
active-region flares. Radio emission associated with XBPs has been obs erved, but in man y 
instances it is consistent with thermal instead of nonthermal emission jNitta et aljfl992h . 
Network flares, associated with the magnetic network boundaries, are about an order of 
magnitude fainter th an XBPs and are more transient, lasting for only about 10 minutes 
dKruckeretaill 1997b . Radio emission has also been observed in these events, but again in 
the majority of the cases it is consistent with thermal emission, and in only a few cases could 
it be associated with nonthermal gyrosynchrotron emission from accelerated electrons. 

Determining whether there is nonthermal emission from a population of accelerated 
electrons outside of active regions and in the quiet Sun would provide important insights into 
the nature of possible small-scale steady-state energization processes in the solar corona. 
Although RHESSI has unprecedented sensitivity over 3 — 25 keV, which has greatly aided 
the study of small active-region flares (see Section l2T2l . emission outside active regions still 
remains elusive. 

Investigating the non-flaring properties of active regions jMcTiernanll2009h or the quiet 
Sun with RHESSI is a non-trivial task. RHESSFs imaging method is designed for bright 
compact flares and ill-suited for the spatially-diffuse emission expected from coronal heat- 
ing. Instead, upper limits to the quiet Sun HXR flux can be obtain by offpointing RHESSI 
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Table 2.2 Table of the power-law indices a found from statistical surveys of the fluxes of X-ray flares. The 
top set are predominantly SXR surveys, the bottom HXR surveys. 



Index a 


Energy Range 


Duration 


Number 


Quantity 


Instrument & Reference 


1.86 


7.7-12.5 keV 


1967(2 weeks) 


177 


Peak 


OSO-3 (Hudson et al. 1969) 


1.75 


1-6.2 keV 1 


1966-1968 


4028 


Peak 


Explorer-33/35 (Drake 1971a) 


1.44 


1-6.2 keV 1 


1966-1968 


4028 


Fluence 


Explorer-33/35 (Drake 1971a) 


1.79 


CA XIX 


1980-1989 




Peak 


BCS (Lee et al. 1995) 


1.86 


1.5-12.4 keV 2 


1980-1989 




Peak 


GOES (Lee et al. 1995) 


1.88 


1.5-12.4 keV 2 


1993-1995 


1054 


Peak 


GOES (Feldman et al. 1997) 


2.11 


1.5-12.4 keV 2 


1976-2000 


49409 


Peak 


GOES (Veronig et al. 2002a) 


2.03 


1.5-12.4 keV 2 


1976-2000 


49409 


Fluence 


GOES (Veronig et al. 2002a) 




1.8 


20keV 


1971-1972 


123 


PeakP h 


OSO-7 (Datlowe et al. 1974) 


2.0 


20keV 


27-06-1980 


25 


PeakP h 


Balloon (Lin et al. 1984) 


1.8 


>30 keV 


1980-1985 


>7000 


Peak 


HXRBS (Dennis 1985) 


1.66-1.75 


>25 keV 


1980-1989 


>7000 


Peak 


BATSE (Schwartz et al. 1992) 


1.61 


>25 keV 


1991 


1262 


Peak 


HXRBS (Schwartz et al. 1992) 


1.75 


>26 keV 


1978-1986 


4356 


Peak 


ICE (Lee et al. 1993) 


1.70-1.86 


>30 keV 


1980-1984 


3578 


Peak 


HXRBS (Bai 1993) 


1.73 


>25 keV 


1980-1989 


7045 


Peak 


HXRBS (Crosbv et al. 1993) 


1.59 


>25 keV 


1980-1989 


2878 


PeakP* 1 


HXRBS (Crosbv et al. 1993) 


1.60-1.74 


>25 keV 


1991-1994 




Peak 


CGRO (Biesecker 1994) 


1.86-2.00 


>26 keV 


1978-1986 


3468 


PeakP h 


ICE (Bromund et al. 1995) 


1.74 


>60 keV 


1980-1989 


12327 


Peak 


HXRBS (Kucera et al. 1997) 


1.58 


>10keV 


1989-1992 


1551 


Peak 


WATCH (Crosbv et al. 1998) 


1.56 


>25 keV 


1991-1994 


5430 


Peak 


BATSE (Aschwanden etal. 1998) 


1.46 


>50 keV 


1991-1994 


5430 


Peak 


BATSE (Aschwanden etal. 1998) 


1.50 


3-6 keV 


2002-2007 


24097 


Peak 


RHESSI (Christe et al. 2008a) 


1.51 


6-12 keV 


2002-2007 


24097 


Peak 


RHESSI (Christe et al. 2008a) 


1.58 


12-25 keV 


2002-2007 


24097 


Peak 


RHESSI (Christe et al. 2008a) 


1.71 


4-8 keV 


2002-2007 


18656 


PeakP h 


RHESSI (Hannah et al. 2008a) 



'2-12A, 2 l-8A. 

The default units for the hard X-ray results are count rate (observed flux) except those indicated by ph which 
are nominally instrument-independent photon fluxes. 



from t he Sun, so that the quiet solar signal is "chopped" as RHESSI rotates llHannah et alJ 
12007 ah . This process has obtained limits that are smaller and cover a wider energy range 
than previously found dHannah et a!]|2OO7bl , l2O10h . as shown in Figure 12.111 A wide range 
of possible thermal and nonthermal emission is still consistent with these limits with only a 
hint as to the relationship between active region and non-active region flares. However, the 
nanoflare model for coronal heating, op erating in a manner similar to active-region flares, 
does seem unlikely dHannah eta"i]|20ld) . Future instruments with higher sensitivity and dy- 
namic range could possibly find such HXR emission, which would be expected from a 
nanofl are heating model even in the absence of flare-like brightenings (e.g., ICargill & Klimchukl 
|2004) . 



3 Flare distributions & parameter scaling 

The occurrence distribution function provides one good way to characterize a large number 
of observations. Many flare distribution functions have been published, and we can sepa- 
rate them into two classes: distributions of raw observable parameters (e.g., peak flux), and 
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distributions of derived parameters (e.g., energy content). A derived parameter such as the 
total event energy might seem a more physically correct object, but the data manipulations 
required to get it will necessarily involve uncertain calibrations and "plausible" assumptions 
about the unobserved parts of the event. To make this latter point more concrete, consider 
the two GOES energy bands (0.5-4 A and 1-8 A). These define an isothermal temperature 
and emission measure, but have little sensitivity at low coronal temperatures (e.g., 1-2 MK). 
Thus in principle a large emission measure and energy can be masked from view and go un- 
accounted for, systematically biasing the energy estimate to be a lower limit. In the following 
we provide a critical review of efforts to understand these different occurrence distribution 
patterns. 



3.1 X-ray flux distributions 

We summarize the reported X-ray flux and fluence (the time-integrated flux) distributions 
in Table l2.2l Most of the power-law indices a fall below the critical value of 2.0, meaning 
that the energetic events dominate, rather than the microflares (information on how to de- 
termine a accu rately is detailed in A ppendix IAY An exception to this is the huge GOES 
event sample of IVeronig et al.l d2002ah . which we showed in Figure [TTT1 but in this case the 
authors deliberately did not attempt to correct for b ackground countin g rates. Indeed, this is 
a fundamentally ambiguous procedure, as noted by |Bornmaniil ( ll990l) . and a typical source 
of systematic error. The SXR observations suffer more from this than the HXR observations 
because of their longer time scales. Determining the fluence also introduces bias: how does 
one determine the time interval for the integration, and how does one allow for sensitivity 
variations (either between instruments or as a function of time)? 

Most of the entries in Table l2.2l (see also the representative plots of several distributions 
in Figure [3~H refer to directly-observed peak fluxes. This is the simplest way to handle the 
data and one that suffers the least from the introduction of systematic errors due to unknown 
or imprecise corrections. Uncertainty due to background rates in the detector still remains, 
as with the intrinsic biases of the sampling. This will have virtually no effect on the largest 
events but dominate over the smallest ones. If the instrument measures spectra in detail, it 
is also possible to take one step away from instrumental bias by fitting a spectrum and then 
evaluating the spectral flux density at a well-measured photon energy, e.g., at 30 keV. 

The representative distributions shown in Figure [3~Tl have several common features: they 
match power laws well over a certain magnitude range, and they roll off towards higher and 
lower magnitudes. The roll-off at the low end is either due entirely to the sensitivity limit 
of the particular instrument, or is heavily confused by selection effects due to this limit. 
In particular, faint events may be missed as their emission is obscured by brighter simul- 
taneous events. A fit to the distribution funct ion that attem pts to characterize this roll-over 
(e.g., a log-normal or a Weibull distribution; |Parnellll2002l) probably has little relevance to 
the physics of the events in the well-observed part of the distribution. At larger magnitudes 
the deficit also might have systematic errors (the saturation of a given detector would be an 
obvious one), but at some point a power-law fit for a <2 will diverge unphysically. Thi s 



implies the existence of a n upper limit of some sort (e.g., iLingenfelter & Hudson! 1980h 



lAlexander & Daoul d2007l) analyzed 10 M- and X-class flares, finding that each was com- 
posed of several bursts that saturated above 20 keV, suggesting a saturation of the HXR 
emission for electron-beams in large flares. 

The sensitivity and bias effects (whether instrumental or observational) result in a clear 
change in the distributions at the extremes, but they can also have a more subtle effect on 
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Fig. 3.1 Several of the flux distribut ions listed in Table l2~2l fcrosbv et alj 1 19931 : [Christe et alj j2 008a: 
ICrosbv et all 1998 ; Hannah e t alj2008al) (left to right, top to bottom). Another example IVeronig et al.l2002at) 
is shown in Figure fTTTI Reproduced by permission of the AAS. 



the power law in the mid-regions. Although the power -law nature of the distribution is not 
dramatically changed, the va lue of Of can be different jAschwanden & Charbonneaull2002l : 
lAschwanden & Parnellll2002h . Some attempts have been made to correct for these biases in 
the context of the derived energy distributions and will be discussed further in Section [3~5l 
The distribution of event energies is obtained by integrating over the spectrum to es- 
timate an energy flux, and then over time. Such conversions are not model-independent. 
Therefore, the direct flux measurement may be a better guide to the general conclusion 
from all such distributions: flare occurrence is scale-invariant. That is, the length scale 
does not change when multiplied by a common factor, a property of power-law distribu- 
tions. Solar flares thus have behavior resembling t hat of earthquakes as described by the 
Gutenberg-Richter law dGutenberg & Richtej|l956h but how general is it for flares? First, 
we note that stellar flare s, on a variety of stellar types, tend to follow similar distributions 
(-e.g. lShakhovs ka]3 ll989h . From region to region, there can be s light variations in t he distri- 
bution, in particular there being a varying upper energy cutoff dKuceraetallll997h . This is 
lower for smaller active regions as there is less free energy, but even in large regions there is 
a finite amount of energy available (making "super flares" not just infrequent but impossi- 
ble). An active region from a period of low solar activity has been found to have a frequency 
distribution which clearly r olls over, this dev iation from a power-law being attributed to the 
low finite energy available i Wheatland 20 1 (t) . 

There can also be variations on intermediate time scales dBail 1 19931 : iBromund et alj 

Il995h . However, in general, these are not strong deviations from the general pattern. From 
RHESSI we find (Figure 13.21 right panel) little evidence for change in the slope of the 
microfiare distribution as a function of phase of t he solar cycle , with e xcellent statistical 
significance. This is consistent with the results of IVeronig et alj d2002ah for GOES events 
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Fig. 3.2 Comparisons of occurrence distributions with phase in the solar cycle. {Left ) Veronig et al. (2002a ) 
using GOES 1-8 A data from solar minimum (1986) to maximum (1989); (right) IChriste et al] I j2008a) using 
RHESSI microfiare data from leaving solar maximum (2002) through to approaching solar minimum (2006). 
Reproduced by permission of the AAS. 



(Figure [3~2l left panel). The scale-invariant property of the flare occurrence distributions is 
thoug ht to provide the evide nce that the coronal magnetic field is in a self-organized critical 
state jLu & Hamiltonin"99 lh and this is detailed further in Section [3~3l 



3.2 Time distributions 

The duration of an event must be known for an estimation of the fiuence and therefore also 
for an estimation of t he event energy. Figu re 13.31 shows two e xamples of duration distri- 
butions: from GOES dVeronig et al.ll2002ah and from RHESSI dChriste et al1l2008ah . The 
GOES data have a clear power-law falloff to long durations (see also lDrakeMl971br) . whereas 
the RHESSI distributions (shown separately for rise, fall, and total times) have more sym- 
metrical distributions. This is because the RHESSI data are considerably more affected by 
selection bias. Trying to determine the duration of the flare crucially depends on being able 
to distinguish the flare from the background. The RHESSI distribution shown is for mi- 
croflares (smaller than GOES C-class) for which it is difficult to separate their start and end 
times from the background. The GOES distribution is only for B-class events and larger 
(>10~ 7 W m~ 2 ) and the background rate in GOES varies less than in RHESSI. 

iLee etal.Ul993h explicitly deal with the simultaneous flux and duration distributi ons of 
an HXR data set, that of (ICE, formerly ISEE-3) above 26 keV dAnderson etaill978l) . This 
analysis of joint variables requires an expl icit consideration of the parameter-space domain 
as a means of understanding sampling bias. lLee ct aD Jl993h concluded that a correlation be- 
tween the duration and peak flux was not present; this would be a requirement if microflares 
or nanofiares were to outweigh the energy present in the ordinary flares. Thus, the power 
laws observed in the distributions of peak fluxes most likely provide a guide to the "true" 
distribution of total event energies. 

The other important temporal signature of flares is the time between one flare and the 
next one, the waiting-time. This is generally taken to be the time between the peak emis- 
sion of one flare and the next. The waiting-time distributions provide information about the 
probability of a flare occurring. 
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Fig. 3.3 Left: the distribution of the durations of flares G OES >B class {Veronig et al. 2002a); right: dis- 
tribution of durations of microflares (Christe et al. 2008a), typically below GOES C class. Reproduced by 
permission of the AAS. 



These distributions can be predicted by self-organized criticality (avalanche) models 
of the coronal magnetic fields, discussed further in Section I3.3I These avalanche models 
have an exponen tial waiting-time distrib ution which corresponds to a constant Poisson flare 
occurrence rate dWheatland et al.ll 1998T) . Such a flare occurrence model produces the de- 
sired power-law distribution of energy and duration. Individual active regions demonstrate 
waiting-time distribution s that are either exponentials or the sum of exponentials dwheatlandi 
l200ll ; lMoonetal.ll200lh . The latter results can be explained by the Poisson occurrence rate 
varying as the active region crosses the solar disk. 

The waiting-time distribution found for large samples of X-ray flares is shown in Fig- 
ure !3.4l The sample of 6,919 flares observed >30 keV with ICE/I SEE 3 is found to have a 
waiting-time distribution that is neither power-law nor exponential dWheatland et al.fl998l) . 
shown left in Figure [341 A larger study of 32,563 GOES C-class and above fla res demon- 
strate s a power-law tail in its waiting-time distribution over long timescales jWheatland 
l2000h. shown right in Figure 13.41 The index of this power law varies with the solar cycle 
( Wheatla nd & Litvine nko l2002h . again consistent with a Poisson occurrence with a time- 
varying rate. An alternative mo del explains this tail using an occurrence rate with a Levy 
distribution dLepreti et al.ll200lh but this also requires a "memory" in the underlying pro- 
cess. This suggests that not only can flare rates be determined but features of the underlying 
physical processes can be understood, although this is still under debate. The waiting-time 
distribution of Type III radio bursts from an a ctive region have been found to be consistent 
with a Poisson process llEastwood et al.ll2010h . 

An overabundance of short waiting-times compared to simulations has now been found 
dWheatland et al.|["l998l) . suggesting that HXR bursts are not inde pendent events. This sym- 
pathe tic flaring behavior has frequently been suggested previously dFritzova-Svestkova et alj 
Il976h . referring to temporally close flares in different active regions. This may also relate to 
the misidentification of several peaks within a single flare as multiple events. It is difficult 
with non-imaging instruments to exclude closely related flares from the same ac tive regions, 
but s tatistically significant evidence for sympathetic flaring has been found dMoon et all 
l2002h . RHESSI would be able to provide such information, as well as investigating waiting- 
times for smaller HXR flares, but the highly discontinuous nature of the data (with gaps 
due to nighttime, South Atlantic Anomaly passage, etc.) would make such analysis highly 
subject to selection effects. 
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Fig. 3.4 Waiting-time distributio ns for different sets of X-ray flares. Left: ICE/ISEE-3 study of 6919 flares 
>30 keV jWheatland et alj|1998l) . the overplotted solid line is a power-law fit. Right: 32,563 GOES flares 
above C-class {Wheatland 2000b, the overplotted lines indicate a model for a time-dependent Poisson process 
(solid) and independent events with an exponential distribution of rates (dashed). Reproduced by permission 
of the AAS. 



3.3 Origin of power-law distribution 

An early analytic model for the statistics of solar-flare occurrence assumed that the avail- 
able energy grew expo nentially with time and would be released as a flare with a Poisson 
probability distribution dRosner & Vaiana|[l978h . This results in a power-law frequency dis- 
tribution. This problem was later recast as a steady-state trans port equation, allowing the 
inclusion of an arbitrary energy resupply rate dLitvinenkdll994h . These models assume that 
all of the free energy is released in a flare, as in a relaxation oscillator. Such behavior has 
never been found in the occurrence patt erns of solar flares , although there are othe r astro- 
phy sical contexts in which it has (e.g., iLewin et al.lll976h . Iwheatland & Glukhovl d 19981) 
and IWheatlandl d2008h generalized these models to take into account that the flares do not 
release all of the free energy by finding a "master equation" to describe the system. This 
model again produces a power-law flare frequency distribution as well as a high-energy 
turnover. 

A complementary view is that the scale-invariant behavior of the flare distributions im- 
plies that the system dynamics can be descr ibed as variations around a self-organized critical 
state dLu & Hamilton! 99 ltlLu et al.l 1993h . In the standard view, magnetic energy builds up 
in complicated and stressed magnetic field structures in the corona due to the motion and 
emergence of magnetic fields through the photosphere. Eventually, the coronal structure 
loses equilibrium catastrophically, and its restructuring (usually thought to involve mag- 
netic reconnection) suddenly liberates some of the built-up magnetic energy. This newly 
released energy goes into the various forms we observe in a solar flare. In the scenario of 
the self-organized critical state, the instabilities are spontaneous, independent of the history 
of how the energy accumulated, and directly trigger a cascade of energy releases, ending in 
a temporarily stable state. Therefore, the same process could easily produce events at the 
nanoflare level as well as the more energetic flares. The avalanche-like behavi or of such a 
system is frequ ently described using an idealised "cellular automaton" model dKatzll 19861 ; 
iBak et al.lfl987h . Such a model does not depend upon the actual physical mechanisms in- 
volved at the microscopic level, replacing them with a schematic set of ad-hoc rules for the 
system evolution. The links between the physics of such a model and its statistical descrip- 
tion remain unclear, though taking a steady-state energ y release in the "master equat ion" is 
similar to assuming an underlying avalanche process I Wheatland & Glukhovlll99"^) . Nev- 
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Fig. 3.5 Distribu tions of HXR power -law index y for Left: 6919 large flares observed above 3()keV 
with ICE/lSEE-3 jBromund et al.lll995t); middle; 55 large flares on the disk (dashed grey) and occulted 
(solid black) with R HESSI iKrucker & Liijl200Sh . Right: 9161 RHESSI microflares from the sample of 
lHannah et al.l l2008al) . the dark red line indicating the error ratio. Reproduced by permission of the AAS. 
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Fig. 3.6 The effect of albedo correction on the observed HXR p ower-law index, with h istograms of y (left) 
found from spectrum fitting and (right) then corrected for albedo I Kasparova et al. 2007). The different lines 
show the power-law fits made over different energy ranges: 398 flares with 70 (solid) 15-20 keV and 7 
(dashed) 20-35 keV, and 123 flares with y, (dotted) 35-50 keV. 



ertheless an extensive lit erature applying cellular automata to solar flares has arisen (e.g., 
ICharbonneau et alj|200ll) . 



3.4 Nonthermal and thermal spectral parameters 

The RHESSI X-ray spectra bridge the thermal and nonthermal spectral domains, at the same 
time providing better spectral resolution than any earlier HXR imager. The nonthermal emis- 
sion of the impulsive phase is often characterized by brief spikes of emission with soft-hard- 
soft (steep to flat to steep) spectral evolution, and its simplest characterization is a power 
law in energy. The SXR thermal emission roughly foll ows the time i ntegral of the HXR 
impulsive-phase emission, following the Neupert effect dNeupertlll968[) . This empirical re- 
sult is thought to demonstrate that the total nonthermal energy deposited heats the chromo- 
spheric plasma, driving it up into flaring loops, producing bright thermal emission. 

3.4.1 Nonthermal 

As regards energetics, the HXR parameter of interest is the slope of the power-law fit to 
the photon spectrum in the impulsive phase, namely 7 in /(e) <x e~T. This spectrum gives 
us the best possible guide to the behavior of the energy in the impulsive-phase electrons, 
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Fig. 3.7 Two examples of low-e nergy break energy di stributions from fitting RHESS I photon flare spectra , 
(left) for 55 large occulted flares iKrucker & Lir]|2008l) and (right) 4236 microflares ftiannah et ai]|2008al) . 
the dark red line indicating the error ratio. Reproduced by permission of the AAS. 

which can be highly significant dBrown|[l97lk iLin & Hudsonll 19761) . Figure [331 shows the 
distributions of the power-law index y from several data sets. Each shows clearly distin- 
guishable selection effects. The left panel of Figure [331 displays the spectral index of large 
flares above 30 keV as found from ICE/ISEE-3 (pre-RHESSI). RHESSI observations of large 
flares indicate similar spectral indices, with occulted flares (those with their footpoints hid- 
den behind the limb s o the HXRs are coronal in origin) have systematically softer spectra 
dKrucker & Linll2008l) as seen in the middle panel of Figure 13.51 RHESSI microflares (Fig- 
ure 13.51 right panel) show many much steeper spectra which are often difficul t to distin- 
guish from the thermal emission (Hann ah et al ] l2008ah . iKasparova et al] d2007l) were able 
to isolate one important source of s ystematic error, namely the albedo resulting from X-rays 
backscattered off the photosphere dSantangelo et alj|l973k iKontar et al.ll2010l) . The albedo 
correction affects mainly low energies, ~20 keV. The result of the correction is a steeper 
photon spectrum at those energies, as shown in detail in Figure [3~6l RHESSI allows imag- 
ing spectroscopy of individual HXR footpoints and the difference in spectral index between 
footpoints in larg e flares (above GOES M -class) with two well-resolved footpoints ranges 
between to 0.6 dSaint-Hilaire et alj|2008l) . 

The essential problem in the interpretation of the hard X -ray spectrum lies in the be- 
havior of the bremsstrahlung cross-section (e.g.. lBrownlll97l|) . The photon spectrum is an 
integral over the electron spectrum at all higher energies. This smearing effect means that, 
in practice, it is difficult to determine the parent electron spectrum at low energies, and yet 
these lower-energy electrons contain most of the energy supplied to the flaring plasma via 
collisions. Analytically, the electron spectrum can be assumed to exist above some "cut- 
off" energy Eq, with means the accelerated electron population appears abruptly at energies 
where the thermal spectrum is negligible. This results in a flattening of the expected photon 
spect ra, with a "break" e nergy occurring at photon energies below the electron cutoff en- 
ergy dKosugi et al]|l988l) . Although the spectral indices o f the source electron distribution s 
and observed photon spectra are related in a simple way (iBrownll 19711 ; IKontar et alj|201(ih . 
this is not the case for the photon break and electron cutoff energies. There are large and 
ill-defined uncertainties in the inverse or forward deconvolutions of the HXR spectrum for 
these parameters. 

ISui et al] d2007t) discuss the behavior of the break energy in a selection of "early impul- 
sive" events, which have less confusion between thermal and nonthermal components. They 
found that 9 out of 33 early impulsive flares demonstrated spectral flattening in the X-ray 
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Fig. 3.8 Histogram of the measured local min- 
imum of the dip in the mean electron flux spec - 
10 12 14 16 18 20 22 trum for a sample of 18 flares fKontar et al.l2008l) . 

Dip Energy (keV) These dips can be removed by correcting for 

albedo and do not represent a low-energy cutoff. 

spectrum towards low energies, a signature of a low-energy cutoff in the electron spectrum. 
In three of these events, albedo correction removed the flattening, but in the remaining six, 
the flattening was found to be consistent with forward-fitting a low-energy cutoff in the 
range 15 — 50 keV, which also correlates with the HXR flux. This is at lower energies than 
found from fitting the X-ray spectrum of flares without early impulsive emission. 

Prior to RHESSI this confusion was not a problem, as an instrumental cutoff of around 
20 — 30 keV meant that the transition of the nonthermal to thermal components of the spectra 
was rarely observed. Figure 13.71 shows two examples of break-energ y distributions derived 
from the RHESSI data, the left panel s howing large occulte d flares dKrucker & Lir]|200o1) 
and the right panel showing microflares dHannah et al.l2008U) . For microflares, the observed 
break tends down to very low energies where there are multiple emission lines in the thermal 
spectra, as can be seen in Figure [2761 This makes it very difficult to determine Eb accurately. 
In addition, a model requiring a sharp cutoff in the nonthermal electron distribution down 
at low energies where the thermal spectrum is non-negligible does not seem appropriate, as 
a smoother transition is expected. The validity of this model and the effort to determine the 
energy in a microflare's nonthermal electrons is problematic and is discussed in detail in 
Sectiorf3~5l 

Instead of trying to fit models to the p hoton spectrum base d upon the expected analytical 
derivation from the electron distribution, iKontar et al.l d200ot) inverted the photon spectrum 
directly to obtain the mean electron spectrum for several flares. They found that there was no 
sharp cutoff in the electron distribution, but instead a dip between 12 — 20 keV. A histogram 
of these dip energies is shown in Figure 13.81 Making albedo corrections to these spectra 
completely removed the dip, suggesting that there is a smooth transition between the thermal 
and nonthermal emissions. Their study also suggests that if low-energy cutoffs exist in the 
mean electron spectra of these flares, they should be located at energies less than ~12 keV. 
This implies that applying a cutoff model is only appropriate when considering an electron 
population well away from the thermal distribution. 

3.4.2 Thermal 

For the basic parameters of the thermal distribution, the present standard approach is to make 
an isothermal fit to a m odel spectrum derived from the Chianti atomic-physics database 
dLandi & Phillipsl2006h . using standard assumptions about abundances and ionization states. 
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Fig. 3.9 Thermal parameter surveys of microflares. (Top left) 291 Yohkoh/SXT microflares i Shimizu 1995) 
(his "active region transient brigh tenings"). Reproduced by permission of the PASJ. (Top right) 9161 RHESSI 
microflares ( Hannah et al, 2008a). The numbered solid lines show expected 4-8 keV RHESSI counting-rate 
levels as a function of model T and EM, the white contour shows the events with GOES T < 10 MK from 
th e bottom panel and th e dashed line (F) the relationship between emission measure and temperature found 
bv lFeldman et alj <1996l) : see Figure l3~10l Note that neither of these studies individually confirms this corre- 
lation, and that the two data sets are almo st disjoint. (Bottom) T he microflare temperature derived using both 
RHESSI and GOES for 6740 microflares ( lHannah et a lj2Q08af). Reproduced by permission of the AAS. 



The isothermal fit determines an effective temperature T and emission measure n 2 V. Fig- 
ure |3.9| shows comparable regions of this parameter space for microflares o bserved with th e 
Yohkoh/SXT grazing-incidence telescope, at energies below about 2 keV (IShimizulll995h . 
top left, and by RHESSI at energies around 6—12 keV jHannah et al .l2008ah . top right. The 
two samples, though taken at different times, arguably represent the same class of events, 
and yet the sets of points are almost disjoint. This illustrates the effects of experimental 
bias, in that the isothermal approximation made for each instrument will produce different 
weightings of the full DEM (differential emission measure) distribution. This can be further 
seen in the bottom panel of Fi gure [3T91 where the temperature has been derived for the same 
microflares at the same time using RHESSI and GOES separately. Clearly each instrument 
is responding to different parts of a DEM distribution. 

Extending this to large flares, the thermal parameters continue to correlate in the sense 
that brighter flares have hig her temperature s , but w ith a relatively slow growth of tempera- 
ture with emission measure. lFeldmanetai] ( ll996l) found it to extend over three decades of 
flare magnitude and a factor of five in temperature (5 to 25 MK) when using GOES emission 
measure and Yohkoh/BC^ peak temperature estimates of 868 A- to X-class flares. But this 
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Fig. 3.10 Thermal fits from various surveys showing the correlation of emission measure (log scale) with 
temperature of solar flare observ ations compared with obs ervations of stellar flares and with quiet-Sun bright- 
enings, from the sources shown (Aschwanden et al. 2008). Reproduced by permission of the AAS. 



correlation cleanly misses the RHESSI microflare points, as shown by the dashed line in 
Figure 13.91 these tend to be much hotter than predicted. The fit misses the SXT parameter 
space in a similar way. We presume these two sampl e s to re present the same physical ob- 
jects, but with differe nt observational bias. iKav et all (120031) found a similar correlation to 
lFeldmanetaT] ( ll996h using GOES emission measure and temperature of 89 B- to X-class 
events. Analysis of super-hot flares (>30 MK) with RHESSI, 37 M- and X-class flares, is 
able to extend parameter space up to 50 MK. The GOES emission of the se events still 
correlates in a similar manner to temperature as in the iFeldman et alJ J 1996h survey, but it 
produces a somewhat flatter correlation in linear-log space dCaspill20ld) ~ 

For individual studies over narrow ranges, there appears to be a poor correlation between 
temperature and emission measure, shown in Figure [3~T0l A possible correla tion appears if 
many studies over a wider range are considered dAschwanden et al.ll2008h incorporating 
even stellar flares and quiet-Sun brightenings. Taken together one seemingly obtains a def- 
inite correlation from sample to sample, if not within a given data set. However, this may 
be force-fitting a single concept to different things; a single power-law fit to all of the points 
would be describable approximately by EM °c T 1 , but if one ignores the quiet-Sun events 
one might prefer a much steeper relationship such as EM k T 15 . Given the strong system- 
atic biases among the different kinds of observation represented, perhaps linked only by the 
word "flare," it is no doubt premature to draw any strong conclusions and caution is required 
when making such comparisons. 



3.4.3 Thermal-nonthermal relationships 

For major flares we note that the impulsive-phase nonthermal signatures often have a simple 
physical relationship with the gradual thermal signatures, namely the Neupert effect where 
the time integral of the HXR (and microwave) emission is empirically seen to match the 
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Fig. 3.11 Comparisons of thermal and nonthermal fluxes for o rdinary flares. Upper: hackground-subtracted 
SXR peak flux from GOES vs. HXR fiuence from HXRBS iVeronig et aI1l2002bl) . The solid lines show 
various functional fits, the dashed lines show proportionality, with extension to smaller fiuence, and the grey 
areas indicate regions below instrumental threshold. The different plot symbols i ndicate the time diffe rence 
between SXR and HXR peaks. Lower: a representative sample of RHES SI fla res I Battaelia et al. 2005), plus 
the same sample but including the enormous stellar flare observed by Osten et al.; (2007J) (the red point). 
Again, the dashed line shows proportionality. 



SXR time profile jNeupedfl968l : lHudsori|[T991al : lDennis & Zarroll 19931) . Does this relation- 
ship extend into the microflare domain as well? Figure [3~TTI shows how well the thermal and 
nonthermal signatures relate. The GOES peak flux and HXRBS (Hard X-Ray Burst Spec- 
trometer on the Solar Maximum Mission, SMM) fluences show a large scatter for weaker 
events, but for the bet ter-observed large en ergetic events approach a proportional relation- 
ship with little scatter dVeronig et aTll2002bl) . One would expect greater scatter in the fainter 
events purely due to selection effects, instrumental sensitivity and analysis procedure (e.g., 

background subtraction performed). 

For RHESSI, a similar effect can be observed llBattaglia et al.ll2005l : lstoiser et al.ll2007t) . 
as shown in the lower panel of Figure [3~TTT l. although here the fluxes are directly compared, 
rather than the SXR flux and HXR fiuence. Also shown is one of the most energetic stellar 
flares observed dOsten et al.ll2007h . over six orders of magnitude brighter in SXRs than the 
largest solar flare. Here again the SXR and HXR fluxes scaled together as in the "big flare 
syndrome." This relationship appears to extend into the microflare domai n. Figure 13.121 
shows that it does, using the large sample of RHESSI microflares studied by lHannah et alj 
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Fig. 3.12 Thermal/nofhermal corre- 
lation for 4236 RHESSI microflares 
I Hannah et al. 2008a). The nonthermal flux 
is the estimated 12 keV emission using the 
fitted nonthermal model. The thermal flux 
is found from fitting the X-ray visibilities 
over 4 — 8 keV. Note that the axes are 
inverted relative to those of Figure l3~TT1 but 
that the line of proportionality is the same. 
Reproduced by permission of the AAS. 



d2008<J) and lChriste etail d2008d) . Note that the deviation from the trend line for the largest 
events is instrumental and due to increased detector deadtime before RHESSI's thin shutter 
deploys. In general, Figure 13.121 provides additional strong evidence that the microflares 
simply r epresent an ex t ension , to lower energies, of the same physical processes at work 
in flares. lAschwandenl d2007l ) has pointed out that such a relationship depends upon the 
scaling of the event environment; at some point we would expect the weakest and the most 
powerful events to show some kind of deviation from this behavior - for example, low- 
altitude loops might have shorter coronal cooling times and therefore depart systematically 
from the observed thermal/nonthermal correlation. 

Detailed quantitative analysis of flare thermal and nonthermal emission suggests that 
the underlying mechanisms are more complicated than empirically shown via the Neu- 
pert effect. The differential emission measures (DEMs) for 80 flares were studied with 
the Yohkoh/SXT and BCS instruments with the finding that the high-temperature plasma 
(> 16.5 MK) is more lik ely to demonstrate the Neupert effect than lower temperature plasma 
I McTiernan et alj[l999h . Given that the Neupert effect is thought to s how that the acceler- 
ated electrons are responsible for heating the chromospheric plasma, IVeronig et al] d2005h 
investigated the similarity of the power in the electron beam compared to the power required 
to produce the observed SXRs. They expected the powers to be better correlated than the 
HXR fluence and SXR flux but found a similar correlation. One possibility for this could be 
that the heating of SXR-emitting plasma is not solely due to the electron beam that produces 
the HXR emission. 



3.5 Energy distributions 

As previously introduced, the thermal and nonthermal energies of flares cannot be deter- 
mined directly from observations but have to be inferred using model assumptions. The 
ambiguities of these models and the errors and bias on the observational parameters result 
in large uncertainties in the energy, and hence in the energy distribution. 
The thermal energy U± can be estimated via 

U th = 3n e k B TV = 3k B T^EM-V ohs f, (3.1) 

where V is the volume of the emitting plasma and / is the filling factor (the ratio of the 
actual volume to the measured value, such that V = /V b s , where n e is the plasma density, 
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T is the temperature, and EM = n 2 V is the emission measure). The temperature and emis- 
sion measure can be found through fitting an isothermal model spectrum to either the flare 
spectrum or the ratio of images in multiple wavelengths. The volume is estimated from the 
2-D area of emission from the images, with some model assumption to convert to a volume. 
The thermal energy (given by Equation 13.1b is an upper limit estimate over the time of the 
obse rved spectrum and images, since a filling factor of / = 1 is normally assumed (see, 
e.g., [F letcher e t al.ll201lh. t hough in reality it could be c onsiderably smaller, 10~ 3 to 10 -2 
dCargill & Klimchuklll997l ; lTakahashi & Watanabell200Cl) . 

The power in the nonthermal electrons of energies above Eq (in keV) is estimated via 

P N (>E C ) = 9.5 xl0 24 r(7- 1)0 (y- 0.5, 1.5)4 _r / ergs" 1 , (3.2) 

where y and Iq are the index and normalization of the power law in the ob served photon 
spectrum (in units of photon spectral flux) and jS(m,«) is the beta function terown|[l97ll ; 
ILir]|l974l) . To estimate the energy in the nonthermal electrons U n the total duration of the 
HXR emission needs to be known. The standard "cold thick-target model" model assumes a 
chromospheric thick target in which a beam of electrons stops, such that the thermal energy 
of the ambient electrons is negligible compared to the energies of the electrons emitting the 
HXR. 

RHESSI crucially allows the thermal and nonthermal energies to be determined in a 
flare using a single instrument, due to the energy range it covers and to its imaging and 
spectroscopic capabilities. Studies of several large flares (C-, M- and X-class) have found 
that the thermal and nonthermal energies are the sa me, to within an order of magnitude 
(e.g.. lHolman et al ] |2003l : ISaint-Hilaire & Bendl2005h . This suggests that the conversion of 
the energy in the accelerated particles into thermal pl asma energy is highly effici ent. For 



C- and M-class events this is typically about 10 3 erg I Saint-Hilaire & Benzll2005h . and at 
minimum about 10 31 erg for large X-class events dHolman et al.l2003L Another study found 
that the power required for the white-light emission was comparable to the electron beam 
power required to pr oduce the HXR emission, but only if the cutoff was less than 25 keV 



power required to pro 
ilFletcheret aDl2007l) . 



For microflares, the situation is more complicated, given the difficulty in determining the 
properties of the photon spectrum and the uncertainty this leads to in estimating the break 
energy. The first detailed analysis on individual RHESSI microflares (A- and B-class flares) 
found that the steepness of their spectra results in considerable power in the th ese low-energy 
electr ons, but only estimated the energy above 25 keV at 10 26 — 10 27 erg jKrucker et al.l 
120021) . In a larger study of microflares, where the observed break found in the photon 
spectrum was used as the cutoff energy, the nonthermal energy ranged over 10 28 ~ 31 erg 
dStoiser et al.l l2007h . This gives an overestimate in the energy since £g< Eq, and using 
slightly larger values for Eq produced typical nonthermal energies of a bout 10 28 erg, similar 
to the thermal energies found. In the largest RHESSI microflare study dHannah et al.l l2008a) 
the thermal and nonthermal energies at about the time of peak emission over 6—12 keV 
were found to be similar, to within an order of magnitude, once an empirical correction 
factor had been used to convert the measured Eb to Eq. This factor was found by fitting a 
broken power laws to model thick-target spectra with a range of indices and low energy cut- 
of fs. Further discussion on the calculation of flare energy and the associated issues is given 
in lHolmaneTal] ( l2010l) . 

The distribution of these RHESSI microflare valu es is shown in Fig ure [2T1 in compari- 
son to previous HXR nonthermal energy distrib utions 1 Crosby et alJl993l : lLin et al. 2001a), 



microflares from the SXR thermal distribution jShimizulll995F and the thermal distribution 
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of EUV nanoflares. dParnell & Juppll200d : lAschwanden et alJl200d : iBenz & Kruckeill2002h 
Although such a figure nicely shows the scaling of nanoflares through microflares to large 
flares, and is crucial for determining coronal heating, it is dangerously deceptive. Each of 
these distributions was found for a different type of event, using various instruments and 
for different periods during the solar cycle, and so each will be affected by d ifferent types 
of selection effect and biases. For instance, the SXT energies dShinuzd[l995h are from 291 
brightenings in o ne active region over five days in 1992 August, whereas the RHESSI mi- 
croflare energies dHannah et alj|20083) are from 9161 events taken over five years, March 
2002 to March 2007. This is a comparison of probably the same type of event but even then 
they are still observed at different wavelengths and over different time periods. Extending 
this comparison to EUV nanoflares is very difficult and may not even be appropriate given 
that we do not know if these are similar events or completely distinct physical processes. 
When looking at the actual values of a that these distributions provide, we see the striking 
observational result that flare p eak flux distribu tions are power laws flatter than a = 2, a 
result that goes back at le ast to | Akabanel Jl956l) . This appears to exclude nanoflares as the 
source of coronal heating dHudsortll991bh . On the other hand, distributions that attempt to 
reflect the eve nt energies themselves have often found steeper spectra that are consistent with 
Parker's idea. lParnell & Juppl l l200Ch found that a ranged from 2.1 to 2.6 f or different model 
assumptions but was consistently in the >2 range. iKrucker&Benzl i ll 9981) found that a was 
between 2.3 and 2.6 but later showed that these values and t hose from other EUV n anoflare 
studies were highly dependent on the flare selection criteria jBenz & Kruckerll2002h . 

So, could selection effects be responsible for this discrepancy between nanoflare and 
larger flare distributions? Attempts h ave been made to recover the intrinsic distributions 
from the observationally derived ones ( lAschwanden et alj200d : lAschwanden & Parnelll2002h 
using a technique similar to that applied to M almquist bias, a threshold selection effect that 
biases galaxy number counts in cosmology dHendrv & Simmon sl ll99d : Iwillickll 19941) . In 
this situation, the counting statistics are biased because the brighter, more distant galaxies 
are more likely to be counted than the fainter ones. This Malmquist-bias procedure has been 
used on the thermal energies derived for EUV nanoflares; it takes the means and covariances 
of the observed data parameters and iterates them back to the intrinsic unbiased values. To 
do this requires model assumptions about how the parameters required to calculate the ther- 
mal energy, in Equation 13.11 relate to each other, so that the probability distribution of the 
intrinsic distribution can be analytically described. However, determining how these param- 
eters scale is, in itself, subject to biases and selection effects, as discussed in Section [3.4.2l 
The analysis procedure may then only adjust the scalings and not test whether the parame- 
ters actually relate in such a manner. Attempting to recover the intrinsic distributions, free of 
instrumental and selection effect biases, should certainly be the priority of any study deriv- 
ing flare distributions, but in practice this does not readily appear to be practicable given that 
the parameter of greatest interest, the energy, is not directly obtainable from the existing ob- 
servations. An alternative approach of determining flare energy, which is much less prone to 
these effects, is to measure the luminous component directly via the resulting change in the 
total solar irradiance (TSI). This has been done for the X17 event SOL2003-10-28T1 1:10 
using the TIM (Total I rradiance Monito r) instrument on the Solar Irradiance and Climate 
Experiment (SORCE) dKopp et al]|2005l) . which is sensitive to the solar emission from X- 
rays to far i nfrared. The estima ted total luminous energy for this flare was approximately 
5 x 10 32 erg dWoods et alj2004T ). This technique is only appropriate at present for the largest 
flares, because of competition from other sources of solar variation. 

Given the evidence so far, it does not appear that flare-like events can heat the quiescent 
corona. Moreover it seems that making the comparison of active-region and quiet-corona 
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flare events is unwise given that th ese transient events hav e different populations in appear 
in different physical environments jBenz & Kruckerll2002h . 

4 Conclusions & discussion 

The RHESSI data have allowed us to study microflares effectively while using the same in- 
strumentation for major events. The statistical study of solar flares has great importance in 
understanding the underlying processes involved in the energy release and subsequent emis- 
sion we observe in the solar corona. Often, emphasis is placed on detailed multi- wavelength 
studies of individual events. Such studies have their merits in revealing insights to the pro- 
cesses in these events. However without being able to place the flare in context of other 
events, it is very difficult to determine whether it is typical or unusual behavior that is being 
studied. Only by studying events in large numbers with a minimum of selection bias can one 
really approach an understanding of the general physics. 

The discovery that "ordinary" flare physics extends down to the tiniest events observable 
by RHESSI (or by GOES) allows us to conclude that such events do not explain coronal heat- 
ing. Even the smallest events RHESSI observes are in the active regions and are flare-like, 
distinguishable as individual events. The nanoflare hypothesis instead requires an apparently 
continuous flare population that is many orders of magnitude smaller. They also continue to 
show the flat peak-flux scaling that puts most of the energy in the most powerful events, 
rather than the weakest ones. In a sense, this conclusion simply confirms the appearance of 
the GOES data - sometimes flare-dominated, and sometimes showing steady emission. This 
is inconsistent with the universality of the flat flare power-law distribution. Of course, nu- 
merous individually unobservable nanoflares could create the apparently steady emissions. 

To make further progress in this field, three crucial things have to happen. Firstly, we 
need spacecraft that have a higher sensitivity, lower background and wider dynamic range 
than RHESSI while maintaining the energy range covered across both imaging and spec- 
troscopy. This would allow the faintest events in active regions and the quiet Sun, as well 
as the fainter emission components of large flares, to be confidently analyzed. One such 
suitable implementation would be a HXR focusing-optics telescope dedicated to solar ob- 
servations. Second, we must understand better how biased our observations are, and how 
we can obtain the intrinsic unbiased physics from our observations. This will require a third 
advance, more sophisticated modeling of these faint events to match improvements in the 
data and their implications. 

Acknowledgements The authors thank the editor (B. R. Dennis) and the two referees for providing many 
constructive comments that have greatly improved this review. IGH would also like to thank the Glasgow SSH 
(Solar Self Help) group for their comments. This work is financially supported through NASA contract NAS 
5-98033, an STFC rolling grant, and by the European Commission through the SOLAIRE Network (MTRN- 
CT-2006-035484), and all are gratefully acknowledged. JK acknowledges support from Grant 205/06/P135 
of the GA CR and the research plan AV0Z10030501. 

A Determining distribution parameters 

For a power-law distribution (Equation ll.U . the standard way to determine the index a is to perform a linear 
fit to the log-log histogram of the data. However this is a highly subjective approach as there is considerable 
choice as to the "best" bin width and fitting method. An alter native and more objective approa ch is to estimate 
the power-law index using the maximum-likelihood method iCrawfordetal.il 9 70: Bai 1993,). This approach 
leads to a remarkably simple calculation on the sample to determine the index above some chosen threshold: 
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N 

a m = —r. h 1. (A.l) 

In this example, the energies U are used, where N is the total number of events and {/, is the energy of the 
i' h event normaliz ed by a threshold energy of Uq. The error in this most likely value of can be estimated 
(Wheatland 2004) as 

a a = (a m -l)N- 1 ' 2 . (A.2) 

The observed distribution, however, is often affected by instrumental and selection effects, resulting in a 
deviation from a power law, so fitting a power law alone would be unwise. A common problem is that the 
smallest events are hard to detect and analyze successfully, resulting in a flattening of the power-law for 
these missing smallest events. To fit this observed bias ed distribution a skew-Laplace distribution can be 
used instead of a single power law {Parnell & Jupp 200(f). This fits the desired power-law distribution for the 
larger events, but below some critical value, a power law with a different index is fitted to the under-reported 
smallest events. The determination of the para meters of this skew-La place distribution again can be found 

using a maximum likelihood estimation method I Parnell & Jupp 2000). 

Another possible distribution to describe the observed sample of data is the Weibull distribution JParnelll 
120021) . This distribution has the form 

/(*K,5)=C(*/$)^ 1) eip[-(*/«) r ]. (A.3) 
where K is the shape parameter and ^ is the scale parameter. For a shape parameter K <1, the resulting 
distribution is similar to a power law, but it turns over at the smallest and largest events. This can then 
represent deficiencies in the smallest and largest events. Again, the under-reporting of the smallest events 
is likely to arise from instrumental sensitivity and selection effect bias. The largest events may be missing 
due to the limited dynamic range of the instrument, if these events saturate the detector; or, there might 
actually be fewer of these events if a critical physical u pper limit is being reached (for instance maximum 
energy available in an active regio n ( Kucera et al. 1997). Such a feature is consistent with the predictions of 
avalanche models iLu et aljl 19931) . as discussed in Section [331 The paramete rs of this Weibull distribution 
can be determined again using the maximum likelihood method )Parnellll2002l) . 

As there are several distributions that could successfully fit the data, a statistical test is needed to deter- 
mine which is best. Such a test is the Kolmogorov-Smirnov statistic (KS), which is the maximu m difference 
betwe en the cumulative distribution function (CDF) and the empirical distribution function (EDF) jPress et a?] 
1992). For a data set of, energies U, the CDF for the i' h energy (/, is the integral of the PDF, using the fitted 
parameters, to (/,. The EDF is derived from the observed/calculated parameters and in this example is the 
number of events with energy less than or equal to Uj, which turns out to be (1 — \/2 )/N. The KS statistic 
then provides a measure of the significance level of each distribution (Press et al. 1992). Plotting CDF versus 
EDF provides a graphical way of determining how consistently the data belong to the chosen distribution. 
A graph similar to the familiar histogram can be obtained by plotting 1— CDF and EDF against the event 
parameters (i.e., energy). 

A detailed example of using the maximum likelihood method to determine the parameters, and testing 
the goo dness of the fit using the KS statistic for power-law and Weibull distributions in the solar context, is 
given in lParnell 42002T) . 
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